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Abstract 

Thermal transport properties of the FPU (3 model with a quadratic pinning term are investi- 
gated for various couplings and temperatures. In particular, the size dependence of the thermal 
conductivity, k oc L", is studied, a agrees with that of the FPU [3 model (with no pinning) at high 
temperatures but decreases at low temperatures. This crossover behavior occurs at a temperature 
depending on the strength of the quadratic pinning. 
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Physics of non-equilibrium is of importance to a broad range of issues in science. Con- 
sequently, non-equilibrium physics has been studied for a long time, yet a number of basic 
problems remain. Thermal transport in one dimensional systems, while seemingly simple, 
has been an area of intriguing active research for a long time[l]. The FPU model is one of 
the most classic models in studies of non-equilibrium physics and a natural model defined 
on a lattice, so that it has been studied from various points of view[l]. By virtue of being 
well studied, many of its properties in the non-equilibrium steady state are known and has 
been quite useful as a guide in studying non-equilibrium physics in various models at finite 
temperatures. One of the key intriguing properties of the model is the lack of bulk behavior 
in thermal transport. Specifically, the thermal conductivity of the model is known to behave 
as K oc L", where L is the size of the system and a ranges from 1/3 to 2/5 l|, |2|, 

In this work, we study thermal properties of a model which generalizes the so called FPU 
/3 model by including a quadratic pinning term and study its properties, from first principles. 
The Hamiltonian of the theory is 



2m 2 ' 2 ' 4 

k=l L -I 

We study the statistical mechanical properties related to thermal transport and analyze how 
the behavior changes with the the physical parameters in the theory. FPU model combined 
with the quartic and quadratic pinning potentials and with randomness have been studied 
recently Isl, flQ. 

Let us explain why we find this model interesting: The FPU model does not have bulk 
behavior, as mentioned above. This can be attributed to the translational symmetry of the 
theory[l|, and is seen in a range of similar models. Any pinning term, such as the one 
we introduce in ([T]) destroys this symmetry, so that the model may achieve bulk behavior. 
However, if we ignore the quartic FPU interaction, the model described by ([1]) is none other 
than a harmonic model whose exact solutions are known and do not display sensible bulk 
behavior {12!. . So by combining these models, we have a situation wherein neither model 



has bulk behavior, yet the combined model might, which we investigate. Also, there are 
arguments as to why the size dependence of the thermal conductivity is universal in theories 
with translational symmetry [2], 0, Q]. As we remove this symmetry, we expect that the size 
dependence will not be the same so that we have a model in which the size dependence can 
change naturally with the parameters of the theory. An objective here is to obtain insight 



into the thermal properties of the model by studying how a changes with the physical 
parameters. Furthermore, the model is a fairly simple one generalizing the well studied 
FPU model. As such, knowing the properties of the model globally in parameter space, we 
hope, will enhance our understanding of the thermal transport properties of similar theories. 

Let us first reduce the model to its simplest form without any loss in generality. Using 
the rescalings in qk^k and time, we obtain the following form of the Hamiltonian, similarly 



to the FPU model 
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Here, the variables are related as q' = u^JrrlfPq, p' = uj'^^/mFJ^ p, jj' = mu'^ij,, 
H' = {rn?'u'^ / I3)H . We also adopted the notation Vg^ = qu+i — qu for convenience. The 
temperature in the two formulations are related as T = /5T'/ (m^cj^), and T is effectively the 
strength of the FPU coupling, similarly to the FPU model and the 0^ theory ^, [l^ . The 
model has a coupling constant, yU, the strength of the quadratic pinning. As we can see from 
the reparametrizations, the quadratic pinning becomes more important for low tempera- 
tures. The size L of the system is a parameter of the model which enters intrinsically into 
transport properties, if the model does not have bulk behavior. We choose to fix the quartic 
coupling and vary T, since it seems more natural to fix the system and vary the temperature, 
but one may equivalently fix T and vary the coupling, due to the above rescaling degrees of 
freedom. 

We now investigate the thermal transport properties of this system by studying its 
non-equilibrium steady states. The energy fiow in the system is defined locally as Jk = 
—Pk^qki^ + Vg^). The system at finite temperature has three parameters fi,L,T and we 
try to elucidate the physics behind the dependence on these parameters, similarly to the 
case of /i = in [3i]. Let us briefiy describe how the thermal conductivity, k, is computed: 
Non-equilibrium steady states are numerically constructed for a given set of /i, L and T using 
thermostats generalizing the Nose-Hoover thermostats, at the boundaries 15|, [iGj. In this 
work, we thermostatted two sites at each end of the system at temperatures (T°, T2*'). Away 
from these thermostatted boundary points, the behavior of the system is governed dynami- 
cally by the Hamiltonian ([2]), including the boundary temperature jumps. One crucial point 
is that the thermostats we employ are known to be able thermalize the harmonic oscilla- 



tor chain 



161 ]. which the classic Nose- Hoover thermostats can not 
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here since as we vary the parameters, the system in some cases will approach the harmonic 
model. The classical equations of motion were integrated numerically using the Runge- 
Kutta method with time steps of 0.002 to 0.02. The physical results were checked to be 
stable against variations in the size of the time step. 10'' to few times 10^ data were averaged 
to obtain the physics results. The system sizes used varied from L = 36 to L = 4000. Fixed 
boundary conditions are used throughout. The standard ideal gas thermometer, T = (p^), 
was used to measure the temperature locally. 

Thermal conductivity is obtained from Fourier's Law 

J = -kVT (3) 

where J is the heat current. To obtain k, multiple non-equilibrium states are constructed 
around the same average temperature T to confirm that Fourier's law is valid and also to 
reduce error. This also allows us to check that the we are not too far away from equilibrium 
and that the linear response still holds. This is performed for fixed values of /i and L. 
Examples of thermal profiles in non-equilibrium steady states are shown in Fig. [H 

In order to obtain k, reliably, the following two points need to be taken care of. Tradi- 
tionally, one often obtains a through the behavior of the heat current J for fixed boundary 
conditions as JL ~ L"^. However, in general, temperature jumps arise at the boundaries, as 
seen in Fig. [1], and this needs to be taken into account [fsl] . Due to the jumps, the gradient 
is not (T9 — T^) /L, but smaller. Since the jumps are proportional to the boundary tempera- 



ture gradients [19|, |20|], they will be smaller for larger systems when the boundary conditions 
are fixed. Consequently, there is a tendency for aj > a to hold. A way to avoid this is 
to measure the thermal gradient inside the system, away from the boundaries, as we shall 
do so here. This is particularly important when the system has parameters such that its is 
close to the harmonic limit and has large boundary jumps. The occurrence of large jumps 
is consistent with the exact analytic results known for the harmonic chains, whose thermal 



profiles consist essentially only of boundary jumps[12|, Il3j. Some profiles with their linear 
fits for deriving the gradient are shown in Fig. [H which illustrate this behavior. Another 
issue that needs to be addressed is the validity of the linear response theory. If we stray too 
far from equilibrium, J will deviate appreciably from linear response theory. Here, we check 
that the linear relation holds and that we keep the gradients to be not too large for each 
parameter set. A priori, there is no rule as to whether J has to be larger or smaller than 
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the linear response prediction, but in the FPU model and the 0^ theory, \J\ is smaller than 



what is expected from /tVT 



2l|. If the boundary temperature differences are fixed while 



varying L, VT is larger for smaller systems so that they are further from equilibrium. This 
effect can lead to an underestimation of n for smaller systems, which in turn can also lead 
an overestimation of a. The problem can be avoided by checking Fourier's law or using 
the knowledge of how large the deviation from the linear response theory is. While it will 
not play an important role here, curvature in the profiles can also affect the extraction of k 
in general. 




FIG. 1: Temperature profiles with /i = 1 for iT^,T^) = (0.8, 1.2), (0.08, 0.12). Temperature has 
been rescaled as 2T/{T^ + T2)- Jumps are larger for the latter boundary conditions. Linear fits to 
the profiles away from the boundaries are also shown but they fit the profile so well that they can 
be seen only close to the boundaries. 



The dependence of the conductivity with various values of fi are shown in Fig. [2l Let 
us try to understand the physics behind the main characteristics of the behavior of k. In 
general, for low and high temperatures, k increases, as in the FPU f3 model (/x = 0). At high 
temperatures, k, for different fi tends to converge since the quartic coupling dominates the 
theory and the effect of the quadratic pinning is relatively small. At low temperatures, the 
theory approaches the harmonic limit and k diverges. Differences due to /i are significant in 
this regime. 

By obtaining k for various values of L with respect to a given set of fi, T, the size 
dependence a may be obtained for the particular set of parameters, fi, T. The dependence 
of a on T for various values of fi are shown in Fig. [31 The behavior of a can be summarized 
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FIG. 2: Temperature dependence of k with L = 132, for /i = (x), /U = 1 (□), fi = i (Q), 
fi = 10 (A) and ^ = 30 (v). The results for the same value of fi have been joined for clarity. 
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FIG. 3: a vs. T for ^ = (x), ^ = 1 (□), /i = 4 (Q), /U = 10 (A) and /i = 30 (v). Inset shows 
the same points (for /i 7^ 0) with the rescaled temperature T/ fi'^ (see text). 

as follows: At higher temperatures, a converges to one value and it is consistent with the 
previous literature for the FPU f3 model jl, 2, sl, Q, 0, [?]. When /i 7^ 0, at lower T, the 
dependence a decreases with T. The crossover behavior to decreasing a occurs at larger T 
for larger fi. This is natural since the effect of the pinning term is greater for larger /i so 
that the crossover occurs at a higher temperature. 

The decrease in a can be understood as the effect of pinning, with the tendency towards 
bulk behavior. This change in a occurs at low T since this is the region where the quadratic 
term is dominant and the quartic FPU coupling is weak. The validity of the argument can 
be confirmed from the dominant term in the potential, which indeed changes as we vary the 



6 



temperature. This crossover behavior can be investigated in equihbrium and an example 
of this is seen in Fig. HI The time averages of the terms in the potential, /ig^, (Vq)l and 
(Vg)^ will equal T, if the term by itself completely dominates the potential, from the virial 
argument. We note here that in the equilibrium system, the expectation values of fiql, (Vg)^ 
and (Vg)^ do not vary inside the system. The above understanding of the crossover allows 
us to compute roughly the crossover temperature, which we denote as Tc. 

r^T {T-> 1), ^iq^ ~ T (T < 1) ^ ~ /x' (4) 

The estimate is consistent with the crossover temperatures seen in Fig. [31 The picture also 
suggests that a more general scaling with the rescaled temperature T//i^ might apply. The 
behavior of a with respect to this rescaled temperature is shown in the inset and the scaling 
works well in practice. We do not have a rigorous argument for this scaling behavior and 
it is worth further study. While there is visibly different behavior below and above this 
"crossover temperature" , we expect this not to be a sharp transition. For T < Tc and a not 
too small, the behavior seems to be reasonably well described by 

T 

a = ao + z/ln— , i/ = 0.055(5) (5) 

J- c 

Here, we denoted the common value of a high temperature constant as oq. If we naively 
extrapolate this behavior to low temperatures, a = is reached at a finite temperature, 
Tcexp(— ao/z/), but we cannot reliably extrapolate to this limit. Paradoxically, numerical 
determination of a is exceedingly difficult for small a values, where we naively expect the 
system to have bulk behavior. This is because in the model described by ([2]), small a values 
occur close to the harmonic limit where the gradients become small and are difficult to 
compute reliably. This difficulty precludes us from definitely predicting what happens in 
the limit a — 0. Several scenarios are possible: The model might reach a = at finite T. 
In this case, if k is finite, the system has bulk behavior. It is also possible that the behavior 
for small T is not described by the log behavior ([5]) and that n diverges in the a ^ limit, 
which is also the T — > limit and this limit coincides with the harmonic model. In this 
case, the model can, in some sense, be arbitrarily close to having bulk behavior yet does not 
achieve it for finite values of the physical parameters. We consider this last possibility to 
be the most likely. It can be noticed that the crossover temperature for the FPU (3 model 
(/i = 0) is Tc = so that it is never reached, if the above scaling is applied. However, the 
model is quite different from 7^ cases so that it should be studied separately. 
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FIG. 4: Temperature dependence of {m^) /T (x), {{Vqf)/T (□) and {{VqY)/T (Q) for ^l=l. 



The behavior of the system as it becomes closer to the harmonic hmit perhaps needs to 
be explained. In the harmonic limit, one might expect k oc L behavior with a being one, so 
that a should rather increase than decrease when T is smaller. This behavior is suggested 
b. .he exaa analv* so.u.ion. fo. the ha™o,.c .„de.0.Q; in these solution, fo. fixed 
boundary thermostat temperatures, T^g, J is independent of L. Naively, the gradient varies 
as (Tg — Tf)/L so that /t should be proportional to L. However, this argument completely 
ignores the boundary effects. With boundary jumps, which are dynamical, the gradient is 
not directly related to the boundary thermostat temperatures and aj = \ > a. In fact, 
in the harmonic case, the gradient is zero inside the system so that k effectively diverges 
even for a finite size system, which is consistent with the behavior of k for small T in our 
results, as in Fig. [2l Given the divergence of k, there is no analytic prediction for a in the 
harmonic limit. One property of the harmonic model is the independence of J with respect 
to L for fixed boundary thermostat temperatures. We find that the dependence of J on L 
does indeed become weaker at lower temperatures. The dynamics of the interior adjusts the 
gradient so that the conductivity has a weaker dependence on L. 

In this work, we computed thermal transport properties of the FPU (3 model with a 
quadratic pinning potential for various values of the pinning potential, temperature and 
system size. We have analyzed the dependence of the conductivity on the physical parame- 
ters of the theory and have obtained an understanding of them. Had we included the quartic 
potential, the theory should have bulk behavior in some parameter region . However, 

in the model discussed, none of the couplings by themselves lead to bulk behavior of the 
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system and the total system does not reach the bulk limit. At low temperatures, where the 
quartic coupling is effectively weak, the quadratic pinning term becomes more important. 
This reduces a (in k oc L") so that the transport behavior is closer to bulk behavior. How- 
ever, the system governed by quadratic terms is harmonic and the conductivity, at the same 
time, approaches divergent behavior. 

The size dependence, a, varies with parameters of the theory. We have used this power 
a in this study but it is worth noting that a priori, there is no rigorous proof that k should 
behave homogeneously as a power of L, even though we find that the behavior applies 
quite well in practice. Some questions remain: In particular, the precise limiting behavior 
of the theory can not be obtained numerically, so that a rigorous theoretical reasoning is 
quite desirable, particularly close to the harmonic limit. To this end, one might envisage an 
application of perturbation theory around the harmonic theory to the FPU quartic coupling, 



similar in spirit to the perturbation theory applied to the quartic pinning coupling [22 1. 
From this perspective, one can view the theory with the quadratic pinning in the FPU 
theory as the extrapolation between the harmonic model with a continuous variation of the 
size dependence of the thermal conductivity. 
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